Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



THE L x - Y x RELATION: USING GALAXY CLUSTER X-RAY LUMINOSITY AS A ROBUST, LOW 

SCATTER MASS PROXY. 

B. J. Maughan 1 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St, Cambridge, MA 02140, USA. 

Draft version February 5, 2008 

ABSTRACT 

We use a sample of 115 galaxy clusters at 0.1 < z < 1.3 observed with Chandra ACIS-I to investigate 
the relation between luminosity and lx (the product of gas mass and temperature). The scatter in 
the relation is dominated by cluster cores, and a tight Lx — Yx relation (11% intrinsic scatter in Lx) 
is recovered if sufficiently large core regions (0.15i?5oo) are excluded. The intrinsic scatter is well de- 
scribed by a lognormal distribution and the relations are consistent for relaxed and disturbed/merging 
clusters. We investigate the Lx — Yx relation in low-quality data (e.g. for clusters detected in X-ray 
survey data) by estimating Lx from soft band count rates, and find that the scatter increases some- 
what to 21%. We confirm the tight correlation between lx and mass and the self-similar evolution 
of that scaling relation out to z = 0.6 for a subset of clusters in our sample with mass estimates from 
the literature. This is used to estimate masses for the entire sample and hence measure the Lx — M 
relation. We find that the scatter in the Lx — M relation is much lower than previous estimates, due 
to the full removal of cluster cores and more robust mass estimates. For high-redshift clusters the 
scatter in the Lx — M relation remains low if cluster cores are not excluded. These results suggest that 
cluster masses can be reliably estimated from simple luminosity measurements in low quality data 
where direct mass estimates, or measurements of Yx are not possible. This has important applications 
in the estimation of cosmological parameters from X-ray cluster surveys. 

Subject headings: cosmology: observations - galaxies: clusters: general - galaxies: high-redshift galax- 
ies: clusters - intergalactic medium - X-rays: galaxies 



1. INTRODUCTION 

The number density of clusters of galaxies and the de- 
tails of their growth from the highest density perturba- 
tions in the early Universe are sensitive to the underly- 
ing cosmology. The high X-ray luminosities of clusters 
makes it relatively easy to detect and study clusters to 
high redshifts with X-ray telescopes. For this reason, 
X-ray cluster studies have been effectively used to im- 
pose tigh t constraints on various cosmolog i cal parame- 
ters (e.g. iHenrvl 120001 IVikhlinin et all 12003 : 1 Allen et all 
l2004f ). Many cluster catalogs have bee n assembled based 
on clusters detect ed by Einstein (e.g. iGioia et al.l|1990f) 
and ROSAT (e g.lScharf et al.lll997l lEbeling et al.l ll998l: 
IVikhlinin et al.lll998HB5hringer et al.ll2001h with further 
surveys plann ed and underway w ith Chandra and XMM- 
Newton (e.g. iRomer et all [2001). These provide large 
statistically complete samples of clusters that are ideal 
for cosmological studies. 

However, the challenge in using clusters as cosmologi- 
cal probes lies in obtaining reliable mass estimates for the 
large numbers of clusters required for statistical studies. 
The most reliable X-ray mass estimates for individual 
clusters are obtained by solving the equation of hydro- 
static equilibrium, which requires measurements of the 
density and temperature gradients of the X-ray emit- 
ting gas. This is possible for nearby, bright clusters, but 
the majority of clusters detected in surveys are too faint 
for such detailed observations. In these cases, cluster 
masses must be estimated from simple properties such 
as X-ray luminosity (Lx) or a single global temperature 
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(kT). Furthermore, statistical samples include many dis- 
turbed clusters for which the assumption of hydrostatic 
equilibrium is undermined, causing additional difficulties 
in deriving their masses. 

Self-similar models, in which the only significant en- 
ergy source in clusters is their gravitational collapse, pre- 
dict simple scaling relations between basic cluster prop- 
erties and the total mass. Observational studies have 
found that such scaling relations exist, but that their 
form is not identical to the self-similar predictions. For 
example, the slope of the X-ray luminosity-temperature 
(Lx — kT) relation is s t eeper than self-similar pr edictions 
[e.g. lMarkevitcblll99a lArnaud fc Evrardlll999f ) and the 
entropy in cluster cores i s higher than predicted (e.g. 
iPonman et"aH [1999, 2003). This indicates the impor- 
tance of non-gravitational effects (such as cooling, merg- 
ers, and jets from AGN) on the energy budget of clusters. 
While the form of the scaling relations differs from self- 
similarity, it is generally found that the evolution of the 
scaling relations follows the self-similar pred ictions (e.g. 
IVikhlinin et~afl l2002t iMaughan etldl I2006D . This sug- 
gests that the evolution is dominated by the changing 
density of the Universe. 

In addition to the shape and evolution of the scaling 
relations, a key consideration is their scatter. Not only 
are low-scatter scaling relations desirable for obtaining 
more precise mass estimates, but the scatter must be 
understood in order to account for the effects of bias in 
samples that are defined based on an observable prop- 
erty (e.g. Lx) that has some finite intrinsic (i.e. not due 
to measurement errors) scatter with mass. The scatter 
in the scaling relations h as been found to be dominated 
by the cluster cores (e.g. lO'Hara et al.ll2006l IChen et alJ 
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2007). The high gas density in many cluster cores leads 
to rapid radiative cooling of the gas which condenses 
and is replaced by gas cooling and flowing in from larger 
radii (Fabian 1994). This runaway cooling is thought to 
be balanced by energy input from AGN or supernovac, 
but leads to bright central peaks in gas density profiles 
and low core temperatures. Such cool-core clusters will 
then deviate from self-similar scaling relations involving 
Lx and/or kT, though these effects can be reduced by 
correcting for the cool-core emission by excluding core 
regions from the measurements (e.g. lMarkevitchlll998T ) 
or including an add i tional parameter to model the effect 
(jO'Hara et al.ll200l iChen et al.ll200l . A second impor- 
tant source of scatter is cluster mer gers, which can caus e 
transient increases in Lx and kT (jRandall et al.l 2002). 
This can cause clusters to scatter along the Lx — kT re- 
lation, b ut away from the mass-temperature (M — kT) 
relation ([Rowley et al J 120041: iKravtsov et al.ll2006l ). A 
further consideration is that the effects of cool cores and 
mergers have different rcdshift dependences. At z > 0.5 
the frequency of cool cores is lower, and o f merging clus- 
ters is higher, than in the local universe (|Jeltema et al.l 
[20051 : iVikhlinin et al.ll2006l : iMaughan et al.ll2007h . ~ 

The X-ray emissivity of cluster gas is proportional 
to the square of its density, so the processes dis- 
cussed above can have strong effects on the luminos- 
ity. Indeed, the luminosity-mass (Lx — M) relation, 
while not commonly stud ied, has been found to hav e 
a large intrinsic scatter ([Reiprich fc Bohringerl l2002f ). 
The M — kT relation, meanwhile has been studied in 
detail and t he scatter is found to be considerably 
smaller Ce.g. iFinoguenov et al.l 120011: ISanderson et al.l 



[200llArnaud et al.ll2005l : IVikhlinin et al . 2006). Finally, 
a cluster's gas mass (M gas ) is believed to be simply re- 
lated to_Jts__totel inass by a universal baryon fraction 
(e.g. I Allen et al.ll2004h . Both M gas and kT are popu- 
larly used as p roxies for the total mass in cosm ologi- 
cal st udies (e.g. IVikhlinin et all 120031: lHenrvi[200l. Re- 
centlv lKravtsov et alj ( 20061 hereafter KVN06) used sim- 
ulations and observations to show that the parameter 
Yx(the product of the X-ray temperature and gas mass) 
is a superior mass proxy to either quantity alone, ex- 
hibiting a low-scatter scaling relation with mass, and be- 
ing insensitive to cluster mergers. This res ult has been 
verifie d independently by the simulations of iPoole et al.l 
(|200l . 

While Yx is an excellent mass proxy, the Lx — M re- 
lation remains of fundamental importance both for esti- 
mating masses when data quality only permit luminos- 
ity measurements, and for characterising the biases in 
cosmological studies based on X-ray flux-limited sam- 
ples. In this paper we use a large sample of clusters 
observed with Chandra to study the Lx — Yx relation 
for the first time, and by using Yx as a mass proxy, we 
investigate the Lx — M relation. A ACDM cosmology 
of H = 70 km s^ 1 Mpc -1 = lOO/i km s" 1 Mpc" 1 , and 
Ojvl = 0.3 (Ha = 0.7) is adopted throughout. All errors 
are quoted at the 68% level. 

2. SAMPLE AND DATA ANALYSIS 

The sample consists of all galaxy clusters observed with 
ACIS-I in the Chandra public archive as of November 
2006 with published redshifts greater than 0.1; 115 clus- 
ters at 0.103 < z < 1.26 (34 at z > 0.5). The con- 



st ruction and anal y sis of the sample is discussed in detail 
in IMaughan et alJ (|2007h . but we repeat some pertinent 
points here. Gas masses were obtained by proj ecting 
the 3D emissivity profile of Vikhlinin et al. (2006, here- 
after V06) along the line of sight, and fitting this to the 
projected emissivity profile computed from the surface 
brightness distribution of each cluster. Temperatures 
were measur ed by fitting an a bsorbed, single tempera- 
ture APEC ([Smith et al.l l2001h model to the spectrum 
extracted from the aperture (0.15 < r < l)i?5oo- Lu- 
minosities were either measured from a spectral fit in 
the aperture of interest, or for the purposes of comput- 
ing values comparable to those measured in survey data, 
from the count rate within that aperture (see ^3.3|) . The 
lower redshift cutoff of the sample ensures that -R500 falls 
within the ACIS-I field of view for all clusters, avoiding 
any extrapolation in the Lx measurements. 

Following the method outlined by KVN06, the radius 
-R500 was determined iteratively, measuring the temper- 
ature in the aperture (0.15 < r < l)i?5oo and the gas 
mass within -R500, computing a new Yx, and hence esti- 
mating a new value of i?5oo- In order to estimate -R500 
fro m Yx we used the Yx — M 50 o relation measured for 
the IVikhlinin et al.l (|2006l ) sample of clusters 
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with B YM = 0.564, C YM = 7.047 x 10 i4 M Q and a YM = 
—2/5 (A. Vikhlinin, priv. comm.). The evolution of the 
relation depends on E(z) = [Om(1 + z) 3 + (1 — Om — 
Sl A )(l + z) 2 + Oa] 1/2 - This observed Y x - M 500 relation 
has a normalisation ~ 15% lower than that found in the 
simulations of KVN06, possibly due to the effect of tur- 
bulent pressure support which is neglected in the mass 
derivations for the observed clusters (KVN06). Adopting 
the normalisation from the simulations would have the 
effect of increasing Lx by a small amount, as -R500 would 
increase by v / 15%. The normalisation of the Lx — M 
relations derived in 551 also would increase by 15%. 

For 4 of the faintest clusters (RXJ0910+5422, 
CLJ1216+2633, CLJ1334+5031 and RXJ1350.0+6007), 
temperatures could not be constrained when the central 
0.15i?5oo was excluded. For those clusters, all of the clus- 
ter emission (r < -R500) was used for the temperature 
measurements in measuring -R500 and Yx- Core regions 
were excluded for the luminosity measurements as re- 
quired. Th e measured properties of all of the clusters 
are given in IMaughan et al 
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3. THE L x - Y x RELATION 

The self-similar Lx — Yx relation can be obtained by 
a simple combination of the M gas — Lx and Lx — kT re- 
lations to give 
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with aLY = 9/5 and Bly = 4/5. This relation was fit 
to the observed Yx and Lx values for the sample, with 
Bly and Cly as free parameters and the luminosities 
divided by E(z) 9 ^ 5 . The best fitting relation was mea- 
sured with an o rthogonal, weighted "BCE S" regression 
(as described bv lAkritas fc BershadvllT996[ ), on the data 



The Lx — Yx relation. 
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Fig. 1. — Lx — Yx relation for bolometric luminosities measured 
in the (0.15 < r < l)iJsoo aperture. Luminosities are scaled by the 
predicted self-similar evolution. 

in log space. The intrinsic scatter was then measured in 
the following way. Consider a set of data (xi,yi) with 
measurement uncertainties (a(x)i,a(y)i) fit by a straight 
line model of slope m and intercept c. The intrinsic scat- 
ter of the data in the y direction (a y ), is defined as the 
value of a y that gives a reduced x 2 statistic of unity. The 
value of x 2 is given by 

With this method, the intrinsic scatter in luminosity, 
a i, was measured for each relation. As the data were fit 
in log space, <jl was multiplied by the natural log of 10 
to give the scatter as a fraction of Lx- This method was 
also used to measure the intrinsic scatter in the x variable 
(Yx, kT or M500). The uncertainty on the measured 
sample was estimated from bootstrap resamples of the 
data. 

The Lx — Yx relations obtained with luminosities 
measured in different apertures are summarised in Table 
[TJ The apertures used correspond to three levels of cool- 
ing core correction, but note that the (0.15 < r < l)i?5oo 
aperture was used to measure the temperature for Yx in 
all cases. The full r < -R500 aperture includes all emis- 
sion from cool cores, and the 70 kpc < r < -R500 aper- 
ture is commonly used to exclude the strongest parts 
of cooling cores, with the resulting luminosity scaled by 
1.06 to account for th e excluded non-cool-core emission 
(e.g. lMarkevitchlll998f l. Finally, the (0.15 < r < l)i? 500 
aperture was chosen to conservatively excluded all cool- 
ing core emission. Bolometric luminosities were used for 
these relations. 

As Table [TJ demonstrates, the scatter in the Lx — Yx 
relation is significantly reduced by excluding cool cores, 
and the (0.15 < r < l)i?5oo aperture is the most effec- 
tive at reducing scatter. The relation derived using this 
aperture is plotted in Fig. [TJ 

3.1. The effect of cooling cores and substructure 

In order to test for any residual effects of cooling core 
emission on the Lx — Yx relation, the sample was split 
into cool core (CC) and non-cool core (NCC) subsam- 
ples. Clusters were classed as CC if the temperature mea- 
sured within r < 0.15i?5oo was cooler than that measured 
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Fig. 2. — Lx — Yx relations for cool core (CC) and non-cool 
core (NCC) clusters. Luminosities are bolometric, measured in the 
(0.15 < r < l)-R5oo aperture, and scaled by the predicted self- 
similar evolution. 

within (0.15 < r < 0.3)i?5oo at a significance of > 2a. 
While this method is dependent on the quality of the 
data used and so is not a robust classification scheme, it 
allows us to separate out the strongest cool core clusters 
to look for departures from the population. The relation 
was fit to the two subsets separately, and the relations 
are plotted in Fig. O There is no evidence for any off- 
set between the CC and NCC populations, though the 
slope of the CC relation is ~ 3cr shallower than the NCC 
relation. The intrinsic scatter of the 18 CC clusters was 
a L = 0.04±0.02, significantly lower than that of the NCC 
clusters (ol = 0.13 ± 0.02). This is likely a reflection 
of the fact that CC clusters are generally more relaxed 
than the rest of the population. When core emission was 
not excluded, the CC clusters were significantly offset to 
higher Lx, as expected. The sample contains clusters 
in a wide variety of dynamical states from relaxed, cool 
core clusters to early and late stage m ergers with a range 
of mass ratios (Maugh an et al.| [2007). To further inves- 
tigate the effect of cluster morphology on the Lx — Yx 
relation, the clusters were split into relaxed and unre- 
laxed subsamples base d on their measured c entroid shift 
{(w)) parameters (see iMaughan et al.ll2007l for details). 
(w) was calibrated by visual inspection of cluster images, 
and the 28 clusters with (w) < 0.005i?5oo were classed as 
relaxed, with the remaining 87 classed as unrelaxed. The 
Lx — Yx relation was fit to these subsets separately and 
the data are plotted in Fig. [3] The normalisations were 
consistent for both subsets, and the slope was slightly 
shallower (~ la) for the relaxed clusters, a similar trend 
to the CC and NCC clusters. The scatter was slightly, 
but not significantly, lower for the relaxed clusters. We 
note that the CC clusters generally have the lowest (w), 
and represent the most relaxed subset of the sample. 

3.2. Evolution 

For the previous fits, we assumed that the Lx — Yx 
relation evolves self-similar ly with a^Y = 9/5. We have 
tested the validity of this assumption by using the data 
to determine the best- fitting value of oly- For different 
values of clly, luminosities were divided by E(z) aLY , the 
best fitting relation was found via BCES regression and 
the x 2 was computed using equation [3] (with a^ = 0). 
The best- fitting qly was the value that minimised x , 
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TABLE 1 

Summary of the best fitting scaling relations between luminosity and Y x , kT, and M500 (the X variables listed in column 1). 
Luminosities were measured in different apertures, either from spectral fits or from count rates as discussed in the text 
and listed in columns 2 and 3. the energy bands used for the soft band luminosities are given in the clusters' rest frames. 
Relations of the form L x = CE(z) a (X/X t ) B were fit to the data with X, = (6 keV,4 x 1O 14 M keV,4 x 10 14 M Q ) and 
a = (9/5, 1,7/3) (self-similar evolution) for lx, KT, and A/500 respectively. Column 4 gives the redshift range of the 

CLUSTERS INCLUDED IN THE FIT. COLUMNS 5 AND 6 GIVE THE NORMALISATION AND SLOPE OF THE L x — X RELATION IN QUESTION, AND 
COLUMNS 7 AND 8 GIVE THE INTRINSIC FRACTIONAL SCATTER IN L x AND X. MASSES WERE ESTIMATED FROM Y x AND THE PARAMETERS OF 
THE BEST-FITTING Lx — M RELATIONS INCLUDE UNCERTAINTIES AND SCATTER IN THE Yx — M500 RELATION. 
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Fig. 3. — L x — Y x relations for relaxed and unrelaxed clusters. 
Luminosities are bolometric, measured in the (0.15 < r < l)-f?soo 
aperture, and scaled by the predicted self-similar evolution. 

and the la confidence intervals were given by the values 
that gave A% 2 = 1 from the minimum. 

The best-fitting evolution was a^y — 2.2 ± 0.1, signif- 
icantly (4<t) higher than the self-similar value. To illus- 
trate the evolution, a local Lx — Yx relation was defined 
by fitting the 20 clusters at 0.1 < z < 0.2. Self-similar 
evolution was assumed for measuring this local relation, 
though the effect at these low redshifts is small. The ra- 
tio of a cluster's measured luminosity to that predicted 
by this local relation for the cluster's Yx is then equal 
to E(z) aLY . Fig. [4] shows the value of this ratio for each 
cluster along with the loci of the self-similar and best- 
fitting evolution. 

While the data prefer a stronger evolution than the 
self-similar prediction, both evolution models give an un- 
acceptable fit (x 2 jv k, 4). This is a result of the signifi- 
cant intrinsic scatter in the data; both evolution models 
require = 0.11 to give a reduced \ 2 of unity. With- 
out a priori knowledge of the details of this scatter that 
would allow us to distinguish between evolution models, 
we find no compelling reason to abandon the theoret- 
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Fig. 4. — Ratio of the observed Lx to that predicted by the local 
L x — Y x relation is plotted against redshift. The curves show the 
predicted self-similar (solid) and best-fitting (dashed) evolution. 

ically motivated self-similar evolution, and continue to 
adopt it throughout this work. 

3.3. Survey-quality data 

The tight correlation between Lx and Yx and that be- 
tween Yx and mass (KVN06) suggests that Lx can pro- 
vide an effective mass proxy. However, data that are of 
sufficient quality (~ 1000 net counts) to allow a spec- 
trum to be fit and a luminosity measured in that way, 
will generally allow Yx to be measured directly (as is the 
case for our sample). The prospect of estimating masses 
reliably from luminosities has much greater potential for 
clusters detected in serendipitous cluster surveys, where 
only count rates are available. We now investigate how 
well the Lx — Yx relation holds up for such relatively 
low-quality data. 

Up to this point, Yx has weakly influenced our Lx 
measurements, as i?soo (and hence the luminosity aper- 
ture) was estimated from Yx- This dependence was 
removed by switching to an aperture of (0.15 < r < 
l)E(z)~ 2 / 3 Mpc, chosen to approximate (0.15 < r < 



The Lx — Yx relation. 
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Fig. 5. — The rest-frame (0.5 — 2.0) keV luminosity measured 
from a full spectral fit is plotted against the luminosity estimated 
from the count rate in the same energy band. The (0.15 < r < 
l)E(z)~ 2 / 3 Mpc aperture was used for both measurements and the 
dashed line indicates a one-to-one correspondence. 

l)i?5oo- Encouragingly, the scatter in the Lx — Yx rela- 
tion was insensitive to this change, with ox = 0.13±0.02. 
Next, the energy band of the luminosities measured from 
the spectral fits was changed from the bolometric band 
to the (0.5 — 2) keV band in each cluster's rest frame. 
This change was motivated by the fact that the soft band 
Lx is less sensitive to the cluster temperature than the 
bolometric band, so can be more reliably estimated from 
a soft band count rate (typical of survey data) without 
a measured temperature. The soft band Lx — Yx rela- 
tion showed slightly increased scatter (ox = 0.21 ± 0.02) 
likely due to the fact that the increased temperature de- 
pendence of the bolometric Lx suppresses the scatter 
somewhat, since Yx is proportional to temperature. 

Finally, for each cluster we measured the net count 
rate in the observed frame (0.5 — 2) keV band, and con- 
verted it to a rest frame (0.5 — 2) keV luminosity. For 
the conversion we assumed an absorbed APEC spectral 
model with the absorption fixed at the galactic value for 
each cluster, the redshift fixed at that of the cluster, and 
metal abundance fixed at 0.3Z Q . The temperature of the 
spectral model was set at an initial guess and Lx was 
calculated. A luminosity-temperature relation was then 
used to estimate kT from L x and the process was iter- 
ated until Lx stabilised. The lMarkevEchl l| 1998ft L x - kT 
relation was used, with self-similar evolution, but the 
conversion is insensitive to the choice of Lx — kT rela- 
tion. The reliability of this conversion was verified by 
comparing our estimated luminosities with those mea- 
sured from the spectral fits, and an excellent agreement 
was found (illustrated in Fig. [5|) . The Lx — Yx relation 
was then measured using these estimated soft band lu- 
minosities, and the scatter did not change significantly 
ox =0.19 ±0.01. 

4. THE Y x - M 50 o RELATION 

The ultimate goal of this study is to investigate the 
Lx — M500 relation using Yx as a mass proxy. The use 
of Yx as a low-scatter m ass proxy has stro ng support 
from simulations (KVN06 IPoole et al.ll2007ft . but there 
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Fig. 6. — Yx — M500 relation using our measured Yx values 
with masses taken from the literature. The masses are scaled by 
the expected self-similar evolution. The line shows the best fit to 
the V06 data; this line is not fit to the data plotted here, but 3 of 
the datapoints are in common with the V06 sample. 

have not yet been many observational studies of this pa- 
rameter. The Yx — M500 relation measured for the V06 
clusters exhibits the same low scatter found in the sim- 
ulated relation (KVN06), but the predicted self-similar 
evolution of the relation has not yet been tested observa- 
tionally. It is therefore desirable to investigate the scat- 
ter and evolution of the Yx — M500 relation using clus- 
ters from our sample. This requires independent mass 
estimates for for the clusters. These can be estimated 
from temperature and gas density profiles measured from 
X-ray observations, assuming hydrostatic equilibrium. 
However, data of sufficient quality are unavailable for 
the majority of our sample, and a full mass analysis of 
the remainder is beyond the scope of the current work. 
Instead, we searched the literature for mass estimates 
for clusters in our sample. We found 12 clusters in our 
sample with values of M500 estimated from full X-ray hy- 
drostatic mass analyses (i.e. using temperature and den- 
sity profiles rather than an isothermal approximation). 
Masses for 3 clusters were taken from the overlap between 
our sample and the V06 sam ple, one was obtained from 
the XMM-Newton sample of lArnaud et al.1 (|2005ft . and 
the remainder came from the XMM-Newton and Chan- 
dra samples presen t ed by fKotov fc Vikhlinin] (|2005f ) and 
iKotov fc~V ikhlinin (|2006ft respectively. The properties 
of the clusters are summarised in Table [2l 

The Yx — M500 relation for these 12 clusters is shown 
in Fig. [6] It was constructed using our measured lx val- 
ues, with the masses taken from the literature and scaled 
by the predicted self-similar evolution (E(z) 2 ^ 5 ). These 
data were compared with the best fit relation to the V06 
clusters (equation [1} , with which three clusters are in 
common. The V06 relation is a good description of the 
data; the reduced x 2 is < 1 so there is no measurable in- 
trinsic scatter. While the intrinsic scatter cannot be well 
constrained with this small sample, the tight correlation 
in Fig. [6] adds further support to the use of lx as a low- 
scatter mass proxy, and provides the first observational 
support for the self-similar evolution of the Yx — M500 
relation to z = 0.6. 

5. ESTIMATING MASSES FROM LUMINOSITIES 

As a final exercise, we use the Yx — M500 relation to es- 
timate masses for the clusters in the sample, and then in- 
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TABLE 2 

Properties of the subset of clusters with masses based on high quality data 
available in the literature. masses were estimated using x-ray hydrostatic 
mass analyses within i?500- the luminosities and y x values presented were 
measured from our analyses, with l x measured in the (0.15 < r < l)i?500 
apertur e. the references for the mass estimates listed in c olumn 6 correspond 

TO: 1 - IVlKHLININ ET AL.1 (120061'! . 2 -IARNAUD ET AL.l (120051). 3 - IKOTOV & VikhlininI 

1120051) . 4 - IKotov fc VikhlininI (120061) . 



Cluster 


z 






M500 


Reference 






(10 44 erg s" 1 ) 


(10 14 M Q keV) 


(1O 14 M ) 




A1413 


0.143 


7.6 + 0.1 


5.6 ± 0.2 


7.8 + 0.8 


1 


A907 


0.153 


4.7 + 0.1 


3.0 + 0.1 


4.7 + 0.4 


1 


A383 


0.187 


3.6 ± 0.1 


1.7 + 0.1 


3.1 + 0.3 


1 


A2204 


0.152 


12.7 + 0.3 


8.9 + 0.7 


8.4 + 0.8 


2 


MS0302. 7+1658 


0.424 


3.0 + 1.2 


9.2 ± 1.5 


2.2 + 0.6 


3 


MS0015. 9+1609 


0.541 


33.9 + 0.6 


13.4 ± 1.0 


8.8 + 1.1 


3 


CLJ1120+4318 


0.600 


8.3 ± 0.7 


2.5 ± 0.6 


4.6 + 1.1 


3 


MACSJ0159. 8-0849 


0.405 


17.7 + 0.3 


11.2 ± 1.1 


9.9 + 3.1 


i 


MACSJ0329. 6-0211 


0.450 


11.3 + 0.6 


3.3 ± 0.4 


3.6 + 1.0 


4 


RXJ1347.5-1145 


0.451 


41.4 + 0.8 


2.4 + 0.2 


13.8 ± 2.8 


4 


MACSJ1621. 3+3810 


0.463 


8.5 + 0.5 


3.8 ± 0.4 


4.4 + 1.0 


4 


MACSJ1423. 8+2404 


0.543 


11.3 + 0.73 


4.1 + 0.7 


4.6 + 1.0 


4 



vestigate the Lx — M500 relation using luminosities mea- 
sured in different ways. In order to include the effects 
of both the intrinsic scatter and the uncertainties on the 
Yx — M500 relation on our derived Lx — M 50 o relation, 
we used a Monte-Carlo approach as follows. A random 
realisation of the Yx — M500 relation (equation [J) was 
created by drawing values of the slope and normalisa- 
tion from the Gaussian distributions defined by Bym = 
0.564 ± 0.009 and C YM = (7.047 ± 0.097) x 1O 14 M . 
These are the best-fitting values from the V06 data, but 
with the errors taken from the fractional errors on the 
KVN06 simulated relation. For each cluster, a mass 
was computed using this relation (assuming self-similar 
evolution) and that mass was then randomised under a 
Gaussian in natural log space centred at the true value 
with a = 0.071. This corresponds to the 7.1% intrin- 
sic scatter in mass in the Yx — M500 relation found by 
KVN06. 

The Lx — M500 relation was then fit using these ran- 
domised masses and our measured luminosities. The pro- 
cess was repeated 1000 times and the distributions of 
measured slope (Blm), normalisation (Clm) and scat- 
ter in Lx — M500 were then used to determine the values 
and uncertainties of those parameters. The uncertain- 
ties and scatter in the Yx — M500 relation did not con- 
tribute strongly to those on the derived Lx — M500 rela- 
tions. In all cases, the ±34 percentiles about the mean of 
the Blm and Clm distributions derived from the Monte- 
Carlo randomisations enclosed a smaller range than the 
statistical uncertainties from the regression fit to the un- 
perturbed Lx — M500 relation. The Monte-Carlo uncer- 
tainties were added in quadrature to the regression uncer- 
tainties to give the final uncertainties for these param- 
eters. We used the median and ±34 percentiles of the 
distribution of <jl from the Monte-Carlo runs as our es- 
timate of the scatter in the Lx — M500 relation. This too 
was not strongly affected by the scatter in Yx — M500; 
the median <jl was <0.03 higher than that measured 
for the unperturbed Lx — M 50 o data in all cases. The 
measured Lx — M500 relations are summarised in Table 
[T] and the relation obtained for bolometric Lx measured 
from spectral fits in the (0.15 < r < l)i?5oo aperture is 




1e+14 1e+15 



M500 (Msol) 

Fig. 7. — Lx — M500 relation for the clusters with masses es- 
timated from the Yx — M500 relation and bolometric luminosities 
measured from spectral fits in the (0.15 < r < l)i?5oo aperture, 
and scaled by the predicted self-similar evolution. 

plotted in Fig. [7] As expected, the scatter follows the 
same trend as in the Lx — Yx relation, reducing from 
0.39±0.04 to 0.17±0.02 when the core regions were fully 
excluded, with a slightly larger scatter of gl = 0.21±0.02 
for survey-quality data. 

5.1. Correlation of scatter in the Lx — Yx and 
Yx — M500 relations. 

Our Monte-Carlo method assumes that the scatter in 
the Lx — Yx and Yx — M 50 o relations is uncorrelated. If 
they were positively or negatively correlated, then the de- 
rived scatter on the Lx — M500 relation would be under- 
or overestimated respectively. Ideally, one would mea- 
sure the offset of clusters in Lx from the Lx — Yx rela- 
tion (Sly) and in Yx from the Yx — M500 relation (8ym) 
and test for correlations between those offsets. This is 
not possible for the full sample, as independent mass es- 
timates are not available, but this can be investigated for 
the subset of clusters with masses taken from the liter- 
ature. The offsets for these clusters are plotted in Fig. 
[HI The data hint at a positive correlation, but the corre- 
lation is not significant if the uncertainties on the data 
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Fig. 8. — Left: Fractional residuals in Lx from the Lx — Yx 
relation (Sly) and m Yx from the Yx — M500 relation (&ym ) for 
the subset of clusters with mass estimates taken from the litera- 
ture. Right: The same residuals for the full sample from one of the 
Monte-Carlo runs for the assumptions of correlated and uncorre- 
cted scatter in the Lx — Yx and Yx — M500 relations. 

points are included. 

In order to investigate the case of positive correlation 
between the scatter in the Lx — Yx and Yx — M 50 o re- 
lations and thus place a conservative upper limit on the 
Lx — M500 scatter, we rederived the Lx — M500 relation 
assuming that Sly and Sym are perfectly correlated. 
In this process, cluster masses were estimated using the 
Yx — M500 relation as before, but were then scattered to 
give a Sym that was proportional to the cluster's known 
Sly- The constant of proportionality used was the ratio 
of the scatter in the Lx — Yx and Yx — M500 relations. 
This is illustrated in Fig. [5] along with an example of 
uncorrelated scatter as assumed previously. With this 
assumption of correlated scatter, the resulting scatter in 
the Lx — M500 relation for Lx measured from spectral 
fits in the (0.15 < r < l)i?5oo aperture increased from 
a L = 0.17 ± 0.02 to a L = 0.23 ± 0.02. For compari- 
son, the Lx — M500 relation was also constructed for the 
subset of 12 clusters with independent mass estimates 
from the literature, using our measured luminosities in 



the (0.15 < r < l)-Rsoo aperture. The scatter of these 
data about the Lx — M 50 o relation derived for the full 
sample was found to be ox = 0.19 ± 0.11. The large 
uncertainties on the scatter due to the small number of 
clusters with independent mass measurements limit our 
ability to pin down the extent of the correlation between 
Sly and 8ym> and the precise value of the scatter in 
the Lx — M500 relation. This will be addressed in future 
work by deriving masses for a larger subset of our sam- 
ple. In the mean time, we continue with the assumption 
of uncorrelated scatter, but the increase in <jl for cor- 
related scatter indicates the size of any sytematic effect 
due to this assumption. 

5.2. The effect of including core emission for high-z 
clusters 

A source of concern for the application of this method 
is that a substantial fraction of the cluster emission is 
excluded by excluding the cluster core out to 0.15i?soo 
(or 0.15E(z)~ 2 / 3 Mpc). To quantify this, we compared 
the flux measured from spectral fits in the (0.15 < r < 
l)-Rsoo aperture and the total flux within -R500 after cor- 
rection for cool cores. The standard cool-core correction 
of multiplying the flux measure d in the 70 kpc < r < 
-R500 aperture by 1.06 was used (|Markevitchl H998L The 
mean ratio of the (0.15 < r < l)i?5oo to total flux was 
0.70, so typically 30% of the flux is excluded. A second 
issue is that for telescopes with a significant point spread 
function (PSF) such as ROSAT and XMM-Newton, the 
central exclusion region becomes smaller than the PSF 
for distant clusters, making core exclusion difficult. 

Recently (jVikhlinin et al.ll2006l ) used measurements of 
the gas density profile slopes in the cores of galaxy clus- 
ters to show that the fraction of clusters with cool cores is 
low at z > 0.5. As cool cores are the dominant source of 
scatter in these scaling relations, it is therefore interest- 
ing to measure the scatter without core exclusion for dis- 
tant clusters to determine if it can be neglected for high-z 
clusters. Luminosities were estimated from count rates 
in the r < lE(z)~ 2 /' i Mpc aperture (i.e. including core 
emission) and the sample was split at z = 0.5 into low- 
and high-redshift subsets. The Lx — M500 relation was 
calculated for each subset as described above, and the pa- 
rameters are given in Table [T] The best- fitting relations 
are consistent for the two subsets, but the scatter for 
the low redshift subset (ax = 0.46 ±0.05) is significantly 
larger than that for the high-z clusters (ax = 0.19±0.05). 
In fact, the scatter for the distant clusters with the cores 
included is the same as that for the full sample with cores 
excluded. This suggests that the core exclusion is gener- 
ally not required for clusters at z > 0.5. Note that while 
the measurement errors are generally larger for the more 
distant clusters, these values for the intrinsic scatter take 
those into account. 

6. DISCUSSION 

Table[T]shows that the scatter in the scaling relations is 
dominated by the cluster core regions. The obvious can- 
didate is the enhanced Lx due to cool cores. However, it 
should be noted that the strongest effects of mergers on 
Lx occur around co re passage of the merging bodies (e.g. 
iRandall et al.l 120021 ) . By excluding a large core region, 
we also effectively remove the strongest part of the scat- 
ter due to mergers. The scatter from the core regions is 



B. J Maughan 



thus due to a combination of mergers and cooling. Once 
the core regions are excised, the morphological status of 
the clusters has little effect on the scatter unless only 
the most relaxed systems (the CC clusters) are consid- 
ered ( §3.ip . The total intrinsic scatter of <tl = 0.36 in 
the Lx — Yx relation can thus be resolved into different 
components. The total scatter reduces by w 70% when 
the central 0.15i?5oo is excluded. A further w 20% of the 
scatter is removed when only the most relaxed CC sys- 
tems are considered; this portion of the scatter is likely 
due to merger effects outside the cores. The remaining 
« 10% of the scatter is likely due to residual substructure 
even in the most relaxed systems, to the scatter in the 
Yx — M500 relation, which has a weak effect on our mea- 
sured luminosities through -R5007 arid to contributions 
from unresolved point sources. 

For comparison with the Lx — Yx relation, we also 
measured the Lx — kT relation for the sample, using lu- 
minosities and temperatures measured in the (0.15 < r < 
l)i?50o aperture. The best fit parameters are given in Ta- 
ble [TJ but the main point to note is that even with our 
most conservative core exclusion, ox in the Lx — kT re- 
lation is three times that in the corresponding Lx — Yx 
relation. The fact that Yx has been shown to be a low- 
scatter mass proxy implies that the larger scatter in the 
Lx — kT relation is due to scatter in the M — kT relation. 
As our mass estimates are based on kT (via Yx), we can- 
not investigate the scatter in the M — kT relation for our 
sample. However, our results suggest that for samples 
such as ours, which include relaxed and disturbed clus- 
ters, luminosity (with a large core region excluded) has 
a tighter relation with mass than temperature does. Al- 
though, recall that Yx has a tighter relation with mass 
than do either Lx or kT. Clusters with strong deviations 
from the Lx — kT relation will be studied in detail in a 
future paper. 

6.1. The shape of the scatter in the Lx — Yx relation 

It is desirable to adopt a functional form for the in- 
trinsic scatter in the scaling relations in order to include 
its effects on e.g. sample completeness and derived mass 
functions. With 115 clusters it is possible to investigate 
the form of the scatter in the Lx — Yx relation. Fig. [9] 
shows a histogram of the log e space residuals from the 
Lx — Yx relation fit to all clusters, using spectral Lx 
measured in the (0.15 < r < l)i?soo aperture. This raw 
scatter is well described by a Gaussian in log space with 
gl = 0.19; a Kolmogorov-Smirnov test gives a null hy- 
pothesis probability of 0.99 that the data are consistent 
with that Gaussian. Here we have implicitly assumed 
that the scatter does not vary with redshift. 

The form of the intrinsic scatter is harder to recover. 
For each cluster we have an observed luminosity (Lx) and 
a luminosity predicted by the Lx — Yx relation for that 
clusters measured Yx {L mo d e i). Each log space residual 
(log e (Lx/ L mo dei)) is associated with a measurement er- 
ror on (Lx/ Lmodei), which we assume is described by 
a Gaussian distribution in linear space. The raw scat- 
ter can then be approximated as a convolution of the 
intrinsic scatter distribution with a linear Gaussian rep- 
resenting the mean measurement error on the residuals. 
This approximation was used to test if the intrinsic scat- 
ter could be described by a lognormal distribution. A 
simulated residual point was drawn from under a Gaus- 




log e (Lx)-log e (L m odel) 

Fig. 9. — Histogram of the residuals from the Lx — Yx relation 
in log e space. This raw scatter is well described by a Gaussian 
of ox = 0.19 (dashed line). The solid line is a Gaussian with 
Cjj = 0.11 corresponding to the measured intrinsic scatter. 
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Fig. 10. — Cumulative probability distributions for the observed 
(solid line) and simulated (dashed line) raw scatter. 

sian with a = 0.11 in log space (our intrinsic scatter 
model). This value of \og e (Lx/ L mo dei) was transformed 
into linear space (Lx/ Lmodei) and then randomised un- 
der a Gaussian of a = 0.17 (the mean measurement er- 
ror on Lx/ Lmodei)- Note that we verified that the mea- 
surement errors on Lx/ L mo dei were not significantly cor- 
related with Lx/ Lmodei- This was repeated for 10,000 
simulated residuals to give a simulated raw scatter distri- 
bution. The resulting distribution was compared to the 
observed raw scatter distribution using the Kolmogorov- 
Smirnov test, and the cumulative probability distribu- 
tions are plotted in Fig. [TO] The distributions agree 
very well, with a null hypothesis probability of 0.91, in- 
dicating that the intrinsic scatter is well described by a 
lognormal distribution. 

6.2. The effect of Eddington bias 

An important factor that we have not yet addressed 
is the issue of Eddington bias. In a flux-limited sample, 
clusters of a range of masses that have fluxes near the de- 
tection limit are scattered into and out of the sample, due 
to the intrinsic scatter in cluster luminosities for a given 
mass. Because the mass function is a decreasing function 
of mass, more clusters will scatter into than out of the 
sample. This results in a bias in both the number of clus- 
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ters and their mean luminosities (as clusters with above 
average luminosities will be over- represented) . The mag- 
nitude of the bias depends on the amount of scatter in the 
Lx — M relation and on the slope of the mass function 
around the flux limit (that is, the mass limit correspond- 
ing to the flux limit when converted to a luminosity at 
the redshift of interest) of the sample. As the mass func- 
tion steepens at higher masses, surveys with higher mass 
limits will suffer a larger bias. The bias can affect the 
slope, normalisation, evolution and scatter of the scaling 
relations measured with the sample. 

The heterogeneous nature of our sample, with clusters 
drawn from a variety of different samples, makes it dif- 
ficult to quantify the effects of the Eddington bias on 
the scaling relations. To illustrate this, consider three 
surveys (all based on ROSAT data) that are represented 
(incompletely) in our sample are the Brightest Cluster 
burvey 7BCS: lEbelm7eT al. 1998) covering z < 0.3, the 
Massive Cluster Survey (MACS: febeling et allboOll ) at 
0-3 < z < 0.6 and t he 400 Square Degree survey (400SD; 
iBurenin et aLll2006f) at 0.3 < z < 1. The flux limits of 
these surveys correspond to masses (within i?5oo) of ap- 
proximately 5 x 1O 14 M (BCS at z = 0.2), 7 x 1O 14 M 
(MACS at z = 0.45), 3 x 10 14 M Q (400SD at z = 0.7). 
The bias thus varies across our sample, and is lowest for 
the deeper surveys. As a simple test of the effect of the 
bias on the measured evolution, the evolution was fit to 
the z > 0.6 clusters alone. These clusters are generally 
drawn from deep surveys for which the bias is lower. The 
best fitting evolution was a^Y = 2.0±0.1, lower than the 
value measured using the whole sample [clly — 2.2±0.1), 
and approaching the self-similar value of 1.8. This indi- 
cates that stronger than self-similar evolution measured 
for the whole sample ( ^3 . 2j) is due in part to the effects of 
bias. However, it should be noted that the evolution is 
measured with respect to a local relation, which is itself 
affected by Eddington bias at some level. The low scat- 
ter in the Lx — Yx relation means that the effects of bias 
should be fairly small, but these will be quantified using 
statistically complete samples at high and low redshifts 
in a future study. 

6.3. Comparison with other measurements of scatter in 
the Lx — M relation. 

The scatter we find in the Lx — M500 relation (ctl = 
0.17 — 0.39) is significantly lower than that measured by 
RB02 who found a total (statistical and intrinsic) scatter 
of 0.74 in a sample of 106 nearby clusters (we refit their 
data to measure the intrinsic scatter alone, and found 
o"l = 0.63). Several important differences exist between 
that work and our study. RB02 used luminosities in- 
cluding the entire core region, derived masses from clus- 
ter temperatures assuming isothermality, and finally, de- 
rived the masses used for the scatter measurement within 
-R20O1 rather than our R§oo- Indeed, for our sample, we 
find ox = 0.39 when all of the core emission is included in 
our Lx measurements (see Table [J), suggesting that the 
remaining difference from the RB02 is due to the scatter 
in the M — kT relation used to derive their masses. 

Other recent work has also shown that the scatter in 
the Lx — M re latio n is lower than pre viously thought. 
iReiprichl (|2006l ) and lStanek et alJ ((2006) compared theo- 
retical mass functions with the X-ray luminosity function 
measured from the RB02 data to determine the scatter 



in the Lx — M relation. The measured scatter depends 
strongly on the assumed cosmology, through the depen- 
dence of the mass function on er 8 and fi m . Using the best- 
fitting cosmological values from the WMAP ye ar 3 results 
(as = 0.74, 9, m = 0.238: ISpergel et al J 12003 1. IReiprichl 
(l2006fl concluded th at the scatter is small (but non-zero). 
IStanek et ail (|2006l l showed that assumin g the WMAP 
year 1 cosmology (erg = 0.9, fi TO = 0.29: ISpergel et all 
2003) results in a large scatter in the Lx — M relation 
and suggested a compromise model with erg = 0.85 and 
Q m = 0.24 that gives a scatter of 0.21 in M200 for a 
given Lx- The measurement from our sample that is 
most directly comparable to this value is the scatter in 
M500 for Lx measured with the cluster c ores included, 
( 7m = 0.21. This is in agreement with the IStanek et al.l 
(2006) "compromise cosmology" value and is wit h in the 
90% upper limit on the scatter set by IReiprichl (|2006l ) 
using the WMAP year 3 cosmology (estimated by con- 
verting their upper limit on the bias factor to a scatter) . 
Note that we do not suggest that the scatter measured 
with our data can be used to distinguish between differ- 
ent cosmological models. 

7. CONCLUSIONS 

We have used a sample of 115 clusters of M500 > 
10 14 M Q spanning the redshift range 0.1 < z < 1.3 to 
study the Lx — Yx relation, and by proxy, the Lx — M500 
relation. We verified the low scatter in the Yx — M500 re- 
lation and its self-similar evolution to z = 0.6 for clusters 
in our sample using mass estimates from the literature. 
Our main conclusions are as follows. 

• A strong correlation exists between Lx and Yx- 

• The scatter in the Lx — Yx (and hence Lx — M500) 
relations is dominated by cluster cores, with a sec- 
ondary contribution from the less relaxed systems 
at larger radii. 

• Once a sufficiently large core region is excluded 
(0.15i?5oo or 0.15£;(z)~ 2 / 3 Mpc) the scatter in the 
Lx — Yx relation, and by proxy Lx — M500 relation 
are low (11% and 17% respectively). 

• The scatter remains reasonably small for survey- 
quality data, where Lx is estimated from soft band 
count rates, increasing to 19% (Lx — Yx) and 21% 
(L x -M500). 

• The slope and normalisation of the Lx — Yx (and 
by proxy Lx — M500) relation are insensitive to the 
merger status of the clusters, although the scat- 
ter can be further reduced by considering only the 
most relaxed clusters. 

• The shape of the intrinsic scatter distribution 
about the Lx — Yx relation is well described by 
a lognormal function. 

• The evolution of the Lx — Yx relation appears con- 
sistent with the self-similar prediction, though our 
lack of knowledge of the details of the scatter and 
Eddington bias limit the strength of this conclu- 
sion. 
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• The inferred scatter in Lx — M500 is much smaller 
than that found by previous studies because of i) 
our improved mass estimates using 13c, and ii) our 
conservative exclusion of core e mission. This is 
con sistent with r e cent results by iReiprichl (|2006l ) 
and lStanek et all (|2006l ). 

• For high-redshift (z > 0.5) clusters, there is no need 
to exclude the core emission. The scatter for high- 
z clusters with cores included is the same as that 
measured for the entire sample with cores excluded. 
This is likely due t o the absence of cool c ore clusters 
at high redshifts (|Vikhlinin et al.ll2006ft . 

These results suggest that a simple luminosity mea- 
surement can provide an effective mass proxy for clus- 



ters with low-quality data (i.e. insufficient to measure 
Yx directly) out to high redshifts. This has important 
applications for testing cosmological models with current 
and future X-ray cluster surveys. 



We thank Christine Jones, Alexey Vikhlinin and 
Thomas Reiprich for useful discussions. We are grate- 
ful to the anonymous referee for their useful comments. 
BJM is supported by NASA through Chandra Postdoc- 
toral Fellowship Award Number PF4-50034 issued by the 
Chandra X-ray Observatory Center, which is operated 
by the Smithsonian Astrophysical Observatory for and 
on behalf of NASA under contract NAS8-03060. 



REFERENCES 



Akritas M. G., Bershady M. A., 1996, ApJ, 470, 706 

Allen S. W., Schmidt R. W., Ebeling H., Fabian A. C, van 

Speybroeck L., 2004, MNRAS, 353, 457 
Arnaud M., Evrard A. E., 1999, MNRAS, 305, 631 
Arnaud M., Pointecouteau E., Pratt G. W., 2005, A&A, 441, 893 
Bohringer H., Schuecker P., Guzzo L., Collins C. A., Voges W., 

Schindler S., Neumann D. M., Cruddace R. G., De Grandi S., 

Chincarini G., Edge A. C, MacGillivray H. T., Shaver P., 2001, 

A&A, 369, 826 

Burcnin R. A., Vikhlinin A., Hornstrup A., Ebeling H., Quintana 
H., Mescheryakov A., 2006, ArXiv Astrophysics e-prints 

Chen Y., Reiprich T. H., Bohringer H., Ikebe Y., Zhang Y. ., 2007, 
ArXiv Astrophysics e-prints 

Ebeling H., Edge A. C, Bohringer H., Allen S. W., Crawford C. S., 
Fabian A. C, Voges W., Huchra J. P., 1998, MNRAS, 301, 881 

Ebeling H., Edge A. C, Henry J. P., 2001, ApJ, 553, 668 

Fabian A. C, 1994, ARA&A, 32, 277 

Finoguenov A., Reiprich T. H., Bohringer H., 2001, A&A, 368, 749 
Gioia I. M., Maccacaro T., Schild R. E., Wolter A., Stocke J. T., 

Morris S. L., Henry J. P., 1990, ApJS, 72, 567 
Henry J. P., 2000, ApJ, 534, 565 
Henry J. P., 2004, ApJ, 609, 603 

Jeltema T. E., Canizares C. R., Bautz M. W., Buote D. A., 2005, 

ApJ, 624, 606 
Kotov O., Vikhlinin A., 2005, ApJ, 633, 781 
Kotov O., Vikhlinin A., 2006, ApJ, 641, 752 
Kravtsov A. V., Vikhlinin A., Nagai D., 2006, ApJ, 650, 128 
Markevitch M., 1998, ApJ, 504, 27 

Maughan B. J., Jones C, Forman W. R., VanSpeybroeck L., 2007, 
ApJ, submitted 

Maughan B. J., Jones C, Jones L. R., Van Speybroeck L., 2007, 
ApJ, 659, 1125 

Maughan B. J., Jones L. R., Ebeling H., Scharf C, 2006, MNRAS, 
365, 509 

O'Hara T. B., Mohr J. J., Bialek J. J., Evrard A. E., 2006, ApJ, 
639, 64 

Ponman T. J., Cannon D. B., Navarro J. F., 1999, Nature, 397, 
135 



Ponman T. J., Sanderson A. J. R., Finoguenov A., 2003, MNRAS, 
343, 331 

Poole G. B., Babul A., McCarthy I. G., Fardal M. A., Bildfell C. J., 

Quinn T., Mahdavi A., 2007, ArXiv Astrophysics e-prints 
Randall S. W., Sarazin C. L., Ricker P. M., 2002, ApJ, 577, 579 
Reiprich T. H., 2006, A&A, 453, L39 
Reiprich T. H., Bohringer H., 2002, ApJ, 567, 716 
Romer A. K., Viana P. T. P., Liddle A. R., Mann R. G., 2001, ApJ, 
547, 594 

Rowley D. R., Thomas P. A., Kay S. T., 2004, MNRAS, 352, 508 
Sanderson A. J. R., Ponman T. J., Finoguenov A., Lloyd-Davies 

E. J., Markevitch M., 2003, MNRAS, 340, 989 
Scharf C, Jones L. R., Ebeling H., Perlman E., Malkan M., Wegner 

G. , 1997, ApJ, 477, 79 

Smith R. K., Brickhouse N. S., Liedahl D. A., Raymond J. C, 

2001, ApJ, 556, L91 
Spergel D. N., Bean R., Dore' O., et. al. N., 2006, ArXiv 

Astrophysics e-prints 
Spergel D. N., Verde L., Peiris H. V., Komatsu E., Nolta M. R., 

Bennett C. L., Halpcrn M., Hinshaw C, Jarosik N., Kogut A., 

Limon M., Meyer S. S., Page L., Tucker G. S., Weiland J. L., 

Wollack E., Wright E. L., 2003, ApJS, 148, 175 
Stanek R., Evrard A. E., Bohringer H., Schuecker P., Nord B., 

2006, ApJ, 648, 956 
Vikhlinin A., Burenin R., Forman W. R., Jones C, Hornstrup A., 

Murray S. S., Quintana H., 2006, ArXiv Astrophysics e-prints 
Vikhlinin A., Kravtsov A., Forman W., Jones C, Markevitch M., 

Murray S. S., Van Speybroeck L., 2006, ApJ, 640, 691 
Vikhlinin A., McNamara B. R., Forman W., Jones C, Quintana 

H. , Hornstrup A., 1998, ApJ, 502, 558 

Vikhlinin A., VanSpeybroeck L., Markevitch M., Forman W. R., 
Grego L., 2002, ApJ, 578, L107 

Vikhlinin A., Voevodkin A., Mullis C. R., VanSpeybroeck L., 
Quintana H., McNamara B. R., Gioia I., Hornstrup A., Henry 
J. P., Forman W. R., Jones C, 2003, ApJ, 590, 15 



